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Abstract 

; A construction of differential constraints compatible with partial differential 
; equations is considered. Certain linear determining equations with parame- 
! ters are used to find such differential constraints. They generalize the classi- 
cal determining equations used in the search for admissible Lie operators. As 
applications of this approach non-linear heat equations and Gibbons-Tsarev's 
; equation are discussed. We introduce the notion of an invariant solution under 
' an involutive distribution and give sufficient conditions for existence of such a 
solution. 
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: 1. Introduction. 

As is well known, one can produce many exact solutions of partial differential 
; equations by means of additional constraints [1], [2]. Differential constraints 
I arisen originally in the theory of partial differential equations of the first order. 
'. Lagrange in particular used differential constraints to find the total integral of 
nonlinear equation 

F{x,y,u,Ux,Uy) = 0. 

Darboux [3] applied differential constraints to integrate the partial differential 
equations of the second order. The detailed description of the Darboux method 
can be found in [1], [4]. 

The general formulation of the method of differential constraints requires 
that the original system of partial differential equations 

F^ = 0,..., F"^ = (LI) 
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be enlarged by appending additional differential equations (differential con- 
straints) 

/ii = 0,..., hp = 0, (1.2) 

such that the over-determined system (1.1), (1.2) satisfies some conditions of 
compatibility. 

The theory of over-determined systems was developed by Delassus, Riquier, 
Cartan, Ritt, Kuranishi, Spencer and others. One can find references in the 
book of Pommaret [5] . Now the applications of over-determined systems include 
such diverse fields as differential geometry, continuum mechanics and nonlinear 
optics. Unfortunately the problem of finding all differential constraints com- 
patible with certain equations can be more complicated than the investigation 
of the original equations. Therefore it is better to content oneself with finding 
constraints in some classes, and these classes must be chosen using additional 
considerations. 

It was recently proposed a new method for finding differential constraints 
which uses linear determining equations. These equations are more general 
than the classical determining equations for Lie generators [6] and depend on 
some parameters. Given an evolution equation 

Ut = F{t,X,U,Ui, . . . ,Un), (1.3) 

where = then according to [7] the linear determining equation corre- 
sponding to (1.3) is of the form 

n i 

Dtih) = E E huD]^\F^^jDr{h), 6,, G R. (1.4) 

Here and throughout Df^ are the operators of total differentiation with re- 
spect to t and X. Equality (1.4) must hold for all solutions of (1.3). The function 
h may depends on The number p is called the order of the 

solution of the equation (1.4). If we have some solution the corresponding 
differential constraint is 

h = 0. (1.5) 

It was also shown in [7] that the equations (1.4) and (1.5) constitute the com- 
patible system. 

The organization of this chapter as follows. In section 2 we introduce the 
requisite concepts and derive a nonlinear determining equation. Linearizing and 
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slightly this equation modifying one we obtain the linear determining equation. 
In section 3 we focus on the solutions of the second and third orders to the 
linear determining equation for Gibbons-Tsarev's equation [8] 

Uu = UxUtx- UtUxx^^- (1-6) 

This gives the corresponding differential constraints and allow us to find some 
exact solutions of (1.6). Section 4 is devoted to systems of reaction-diffusion 
equations 

Vt = di{u'^v''vx)x + f2{%v). 

In the final section the invariant solutions under an involutive distribution are 
discussed. We consider the problem of finding involutive distributions that 
enable us to obtain invariant solutions to evolution equations. 



2. Linear determining equations for differential constraints. 

We consider a system S of evolution equations: 

ui = F\t, x,u\..., vT^, ul.ul,...), i = 1, . . . , m. (2.1) 

The right-hand sides of (2.1) can depend on the derivatives of arbitrary orders 
with respect to x. For simplicity we assume that the F* are smooth functions. 
We supplement S by differential constraints H: 

hj{t,x,u\...,u'^,ul,ul,...) =0, j = l,...,p (2.2) 

where p < m. The functions hi can depends on derivatives of arbitrary orders 
with respect to x, but contain no derivatives. 

According to [11], the differential constraints H define an invariant manifold 
of the system S if 

Dtihj)iis]n[H] =0, j = l,...,p. (2.3) 

We denote by [S] the equations (2.1) and their differential consequences with 
respect to x. The constraint (2.2) and its differential consequences with respect 
to X are denoted by [H] . 
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If there are m differential constraints which are all resolved with respect 
to highest derivatives of all functions u^,...,u^, then, as shown in [18] the 
invariance condition (2.3) yields the existence of a smooth solution of the system 
(2.1), (2.2). 

For simplicity we now consider an equation of the second order 

ut = F{t,x,u,ui,U2) (2.4) 

and the differential constrain 

h{t,X,U,Ui, ...,Un)=0, 

where Uk = We denote by [E] the equation (2.4) and its differential conse- 
quences. 

The invariance condition (2.3) is equivalent to the following nonlinear equa- 
tion 



+F,Mu^-.u.., - 2FMhu^.j)h, (2.5) 



with n > 4. 

Indeed, it is easy to see that 



A(M|[i.] ^ D-,iF) + K^_,D--\F) + K^_,D--\F). (2.6) 

Here and throughout we write a ~ to indicate that there are no terms 
including Wn, Wn+i, Un+2 in the difference a — j3 . Since n > 4, the terms on the 
right-hand side of (2.6) can be represented as follows: 



Dl{F) ~ Fu,Un+2^Fu,^nD,{Fu,)]Un+l+ 



Un{Fu^ + (n - l)Dx{FuS) + 



h D'^^'^(F) II F h 

J — '^n-^ U2"'Un-2- 
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Hence (2.6) can be written as the relation 



Dt{h)\[E] ^ Fu^Un+2 + Un+i[Fu^ + nDx{Fu^) + Fu^hu^_^] + 



Fu+ 



n{n — 1) 



Un-2 



It is easy to see that 

Dl{h) ~ + Un+lhu^_, + + 2Dj;{hu„_J - Unhu^_^u^_^]. 

Hence the difference 

Dt{h\E] - FuM{h) - [Fu,+nD,{Fu,)]D,{h) 

contains no terms with Un+2 and Un+i- A direct calculation shows that there 
are no terms containing Un in the expression for the function 



where 



l = Dt{h\E]-M{h)., 



M{h) = Fu^Dlih) + [Fu, + nD,[Fu,)\D,{h) + 



Fu + nD,{Fu,)- 



-K^_,D,{Fu,) + 



that is, 7 ~ 0. We claim that 7 is equal to zero. Since H is an invariant 
manifold, it follows that 

M{h) +7 = 

on the set [E] n [H]. Clearly, M{h) vanish there, therefore 

i\[EnH] = 0. 

Since 7 is independent of u^, w^, . . ., we can rewrite the last equality as follows: 

7|[i?] = 0. 

As shown above, 7 ~ 0, that is 7 can depend only on Un-i, Wn-2, • • • • On the 
other hand, h depends on Un- Hence 7 is equal to zero. 



It is clear that the equation (2.5) is difficult to solve, therefore, in place of 
the non-linear equation, we propose to use in the search for invariant manifolds 
linear equations of a similar kind. This leads to the following definition. 

Definition. The equation of the form 



[c^Fu + c^D,{Fu,) + c^DiiFjjh (2.7) 

is called the linear determining equation corresponding to (2.4). Here Ci, . . . , C5 
are some constants. 

It is easy to generalize the previous definition for case of equation (1.3). 

Definition. A linear determining equation corresponding to (1.3) is the equa- 
tion (1.4). 

The nonlinear diffusion equation 

Ut = {u^Ux)x + f{u) 

have been considered in [12]. It is shown that the solutions of linear determining 
equations of the second and third orders exist only if / belongs to the special 
forms. 

In particular the equation 

Ut = {u~^^'^Ux)x + 'ITT'U — 2ky/u, 171, k G R 

is compatible with differential constrain 

- ^ + ^ + re-'^'/\'/' + se^'/'V^ = 0, r,seR. 
2u Aw^ 



This leads to the following representation 

2X'T' 



u = 



(X + T)2 

The functions T{t) and X{x) satisfy the ordinary differential equations 

{X'f = (csX^ + C2X^ + ciX + co)^ 
(T'f = Ai-csT^ + C2T^ - ciX + Co)', 



where 03,02, ci and cq are arbitrary constants, A = (— 2r)^/^. It is possible to 
show that the functions X and T can be expressed by means of elliptic functions 
[12]. The others examples can be found in [12]. 

The linear determining equation corresponding to the equation 

Utt = F{t,X,U,Ui, . . . ,Un) (2.8) 

is as follows: 



Tl Z 

Dtm^E] = E E bikDi-\F^^_,)D:-\h), hu e (2.9) 



i=0 k=0 



where [E] means the equation (2.8) and its differential consequences. 
Example. Let us consider the equation 

ah . / X 

utt = Uxx + T^t + -'^x + sm{u) , (2. 10) 

t X 

It is easy to check that the symmetry group [6] for (2.1) is trivial. On the other 
hand, the linear determining equation 

Duh = D^^h + ^Dth + -D:,h + (cos(w) = 

t X X'^ 

corresponding to (2.10) has the solution h = xut + tUx- Differential constraint 

h = xut + tUx = 

gives the following representation 

u = V{x'^-t^) (2.11) 

for a solution of (2.10). Substituting (2.11) into (2.10) we have the ordinary 
differential equation 

AzV" + 2(6 - a + 2)V' + sin(y) = 0. 



3. Gibbons-Tsarev's equation. 
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In this section we will consider the Gibbons-Tsarev equation [8] 

■^xx ~l~ ^y^xy ^x^yy ~l~ l 0- (^■"^) 

which arises in reductions of the Benney equation. 

By (2.9), the linear determining equation has the form 

D^h + ZyDxDyh — ZxDyh + hiZyyDxh + b2ZxyDyh = 0, (3.2) 

The constants 61 and 62 are to be determined together with the function h. 
It can be shown that the equation (3.2) has a solution of the form 

h = Zyy + Q(^X, y, Z, Zxi Zy, Zxx) 

if and only if the function g is independent of Zxx- Therefore we shall start with 
solutions of the second order 

h = Zyy + g{x,y,z,Zx,Zy). (3.3) 

Substituting (3.3) into (3.2) leads to an equation which includes derivatives of 
the third order. We can express all mixed derivatives by means of (3.1). Setting 
the coefficients of Zxxx and Zyyy equal to zero we obtain : 

61 = 1, 62 = -1. 

The left-hand side of (3.2) is a polynomial with respect to Zxx and Zyy. This 
polynomial must identically vanish. Collecting similar terms we have the fol- 
lowing equations: 

P9pp + qgpq - 9qq + 2gp = 0, (3.4) 
{-q^ - 2p)gpp + 2gqq + q^gxp + q{q^ + 2p)gyp + q'^{q'^ + Sp)gzp- 

-2qgxq - q^Qyq " Q^Q^ + 2p)^^g + qgy + q^gz - ^Qp = 0, (3.5) 
2p^gpp - {q^ + 2p)gqq + pq'^gxp - 2p^qgyp - p^q^gzp+ 
+q{q^ + 2p)gxq - pq^gyq + 2p^qgzq + q^gx - pqQy + ^PQp = 0, (3.6) 
P9pp - 9pq - 9qq + (1^9xp " 2pqgyp - pq^gzp + 2qgxq + q^gyq + q{q^ + 2p)gzq- 

-q^{q^ + 2p)gxz+pq^gyz-p^q^9zz + 2gp-q^gxx-q^9xy+pq^9yy-(19y = 0, (3.7) 
where p = Zx and q = Zy. 



It is possible to show that the general solution of the equations (3.4) - (3.7) 

is 

h = Zyy + ci{z^ + {Sza: + 4a;) 2:^ + 3yZy + {z^: + 2xf + 2z)+ 

+C2{Zy + {2zx + 3x)zy + 2y) + cs{zy + Zx + 2x) + C/^Zy + C5. 

Hence the differential constraint /z, = is compatible with the Gibbons-Tsarev 
equation (3.1). In the case ci = C2 = C3 = we obtain the differential constraint 

Zyy + C4Zy + C5 = . (3.8) 

From (3.8) we find the following representation 

z = ai exp(-C4?/) - C5?//c4 + 02, (3.9) 

where ai and 02 depend on x. Substituting (3.9) into (3.1) we derive two 
ordinary differential equations 

a2 + 1 = U, fli + c^tti — c^aia2 = U. 

The first equation has the solution 

a2 = —X /2 + cqx + C7. cq, C7 e R 

In this case the second equation is 

a'l + c^a'i + cl{x — cojai = 0. 

Setting ai = exp{—C5x/2)v{x) we obtain the equation 

v" + {A + Bx)v = 0, A,B eR. 

According to [9] the solutions of last equation can be expressed in terms of the 
Airy functions. 

It can be shown that the linear determining equation (3.2) has the following 
solution of the third order 

h = Zyyy + Cii^z^ + (lO^;^; + 12^)4 + Qyzl + {Qzl + ISxZa: + 2z + 12x^)zy+ 

+4yz^ + 6xy) + €2(5^ + {I2z^ + 15x)zJ + 6yZy + 3zl + lOxz^, + 15/2x2+ 
-\-z) + C3{2Zy + {3zx + 4:x)zy + ?/) + C4{3Zy + 2zx + 3a;) + c^Zy + cq. 
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The corresponding constants bi and 62 in (3.2) are given by 

61 = 2, 62 = -2. 

In the case ci = C2 = cs = cq = and C5 = — 1 the function /i gives the 
differential constraint 

^yyy — Zy = 0. (3.10) 
From (3.10) we obtain the following representation 

z = si{x) + S2{x)ey + S3{x)e~y. 

The functions si(x), S2(x) and S3 must satisfy the equations 

§2 - s[s2 = 0, 

Si — 2S3S2 — 2S2S3 + 1 = 0, 

S3 - s[s3 = 0. 

If S3 = as2 then the last system reduces to the two equations 

II I III . 

S2 — S]^S2 = 0, s^ — 4as2S2 + 1 = 0, a E R. (3.11) 

Integrating the second equation, we find that 

s[ = -t -b + 2asl, b e R. 

We can insert this expression in (3.11) and obtain the second-order equation 

S2 + (t + 5 - 2as2)s2 = 

Using the transformations ti = t-\-b and w = y/as2^ we take the equation in S2 
to the second Painleve equation [10] 

w = 2w + tiw. 

The differential constraint 

^yyy ~ ^ 

leads to a solution of (3.1) which is expressed in terms of elementary functions. 



4. The reaction-diffusion equations. 
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In this section we consider systems of the second order 

ut = iu''v^Ua:)x + fh (4.1) 

Vt = di{uVv,),^h, (4.2) 

where I, n, and di are arbitrary constants, the functions /i and /2 depend 
on u and v. This model plays important role in various applications. 

Using the method described in the section 2 it is possible to obtain the linear 
determining equations to the system (4.1) - (4.2). We give the equations in the 
final shape: 

Dt{h) = u^v^Dl{h) + (biku'^-Wu, + b2M-^V:,)D,{h)+ 

+hlu^V^-\^D:,{fi) + (64/1^ + (64 + 265 + hQ)k{vl'-\^U^:,+ 

+ {k - l)u^~'^v^ul) + (64 + 365 + hQ)klv!'~^v^'^u^Vx + (65+ 
+h)l{vl'v'-^v,, + {1- l)u'v''\l))h + (67/1.+ 

+b8liu''v^~'^u^x + ku'^-W-^ul + (/ - l)u''v^-\^V:,)^(3, (4.3) 

Dt{(5) = diu'^v^'Dlif]) + di{bQmu'^-'^v''u^ + bionu'^v''-'^v^)D^{(]) + 
+biidimvJ^-\'^v^D:,{h) + (612/2U + bi^dim{u'^-^v'^V:,:,+ 
+nu'^~^v''~^vl + (m - l)u'^~'^v'^UxVx)^h + (614/2^; + din{bu-^ 

+26i5 + biQ){yJ^v''-\,, + (n - l)vJ^v''-\l) + c?i(6i4 + 36i5+ 

+biQ)'mnu'^^^v^~^UxVx + dim{bi5 + 6i6)(ii'"""^t'"t'Ua;a;+ 

+(m - l)w"^- Vw^))/?. (4.4) 

Here h and are required functions, bij are some constants. 

At first we seek for solutions of the second order to (4.3)- (4.4) of the form 

h = Uxx + hi{t, X, U, V, Ux, Vx), P = Vxx + A(^, X, U, V, Ux, Vx). 

We do not itemize the systems (4.1)- (4. 2) which lead to the simplest solutions 

h = Uxxj P — "^xx 
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of the equations (4.3)-(4.4). Moreover, we did not include solutions generated 
by the symmetry groups. 

The full list of nonlinear systems and corresponding solutions of linear de- 
termining equations is as follows: 



Ut = iuUx)^ + 2/iiW + fuV + /i3, hi = Uj::, + /n, 

Vt = diVxx + J21U V + /23, Pl=Vxx + -r'^ 

Jll Jl2 



Ut = iuUx)x + 2/iiU + fuV + /i3, h2 = Uxx + /ll, 



vt = di (uvx) X + /21W + {Sdifn - 2/12)1; + /23, P2 = v. 



XX 



/12' 



/ N , /ll 2 , ^ /21/12 , , fll 

^ Jll <J 

Vt = -2{uVx)^ + f2lU-fnUV--j-U^, /^3='J^a;ar + ^-^ /ll/l2^ + ^w2- 

^/l2 ^ 0/12 \ 



-/21/: 
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Ut={uU:,)j; + fnUV-\-fi2u'^-\-fl3U + fu: = Uxx + + + 

2/12 Odd 

i;^ =-2(ui;a;)a;- -— (/iiwi;+/i2w2+ ^4 = 0; 

Jll 

+/13W + /14), 

= + /llW + /I5 = + /n, 

Vt = {hlU + h2-^Jl{v))/f[{vl ^5=^,,+ |g^2 + ||-; 
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Ut = {uUx)x + hlU + h2V + /l3, ^6 = Uxx + /ll, 



/l2 /l2 

+/22?^ + /23, 

, r /22/1I , . , , 2/22 

Ut = Uxx + JnU 7 V + /i3, = Uxx + 

/2I <J 

2/21 

Vt = di{uVx)x + dif2lU + dif22V + f23, l37 = Va ' 



3 ' 



Vt = di [vVx)x + /21W + /22^^ + /23, /^S = Vxx " 



= + /uW + /i2^; + /l3, /iS = Wa;ar + ^ .2 ^ (/2]/l2 " /22/11) , 

/21/12 — /22/1I 



3/iic?i 

/ X , r , r , f , , /12/21 - /11/22 

= [VUx)x + /UW + }l2V + /l3, /l9 = Uxx H ; 

<J/21 

, ^ /12/2I - /11/22 

Vt = diVxx + f2lU + f22V + f23, P9 = Uxx ; 

'j/22 

3(c/i/ii-/2i) ' 
/12/21 — /11/22 



Ut = {uUx)x + hlU + fl2U + /l3, hiQ=Uxx 
Vt = di {uv^)x + /2lli + /22^^ + /23, ^10 = U. 



XX 



3(6^1/12-/22) ' 

Ut = ivux)x + fnu + /12?; + /i3, hii = Uxx--r-r 

"1/11 — <J 722 

, / N 2/22((il/ll - 3/22) ^ . f22 . r . r 

Vt = di{vVx)x ^5 Pll = "^^a^x + -1- + h2V + /23; 

"1/12 «1 

/ \ , r , 5/11/22 , f 7 , /11/22 

= (i^W^jx + /llW H 7 V + hi2 = Uxx+^ , 

hi J. ^J21 

Vt = /21W + /22^^ + /23, P12 = Vxx 
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/ I \ , r . 5/11/23 I . r 1 . /11/23 



Ut = (UU^)^ + /uW + /l2f + /l3, /il4 = Uxx: 

Ut = iuUa:)x + /ll^ + /i2?^^/^ + /l3, ^^15 = U^x, 

= {U^^^U^):, + jwU^I^V + /i2W^/^ + /13W, /ll6 = M^x " "T^, 

4u 



= (il(^^x)x + /21^ + /22W^/^ + /23, Ae = 



XXI 

Ut = {u'^/'^U:c)x + /llW-^/2^ + /l2W-^/2 + /13U, /ii7 = ""^^ + 

i;^ = di{vv^)^ + /21?; + /22^i^/^ + /23, A? = ^x^; 

Ut = Uxx + /llW + /l2^^^/^ + /l3, his = 

1;, = di{v'/\,), + /2ini;3/4 + ^^^^3/4 + ^^^^^ = - ^; 

Ut = Uxx + /llW + h2V^''^ + /l3, /il9 = Uxx, 

= di{v^l'^Vx)x + /2lW^;-^/2 ^ /22^-l/2 + J23^^, = ^xa: + ^] 



Ut = {u^l'^Ux)x + hlU^'^V + /i2w3/4 ^ j^^^^ /i20 = W 



20 — U^xx — 



4u ' 

= diVxx + /21^ + f22U^'^ + /23, /?20 = Vxx] 

U^ 

Ut = iu^/'^Ux)x + hlU-^l'^V + fviU-^l'^ + /13W, /i21 = ""^^ + 
= diVxx + /21^^ + /22W^/^ + /23, /?21 = ^^xxi 
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Vt = diiv'N,), + /2inV4^3/4 + ^^^^3/4 + ^^^^^ ^ _ 



vt 



= diiv'/'^V,), + f2lU^I''v^l^ + /22^'/' + /23^, ^^24 = ' 



3i; 



4^; 

^2 

Ut = (wVV)a; + /llW-V2|;l/4 ^ /^^W'^/^ + /13W, /i25 = U^^x + — 



= (uV^w^)^ + fnu^/^v + /i2u3/^ + fi3% h26 = u^x - -r^, 

= + f2lV + /22W^/^ + /23, /?26 = ^a;a;; 

= {v'^/'^U:,)x + /llU + /l2^;^/^ + /l3, /i27 = 

S7;2 

3^2 

Ui = + fnU^^'^V + /i2W^/4 + /13U, /i28 = " "T^, 

4u 

i;^ = c^i (u V V ) . + /2i^ + /22W + , = ; 

Ut = iv^^'^Ux)x + /llW + fuV^/'^ + /l3, /i29 = Uxx: 

3^2 

^2 

Ut = iu^^'^Ux)x + /llW-l/^U + /l2W-l/^ + /13U, /l30 = Uxx + ^, 

= di{u^/h^)^ + /2li; + /22^i^/^ + /23, /?30 = ^xxi 
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Ut = iv^/'^Ua:)x + fllU + fuV^^^ + /13, 



^31 — Uxx: 



Vt = di{v^/^Vx)x + hlUV-^''^ + f22V-^'^ + /23^, /?31 = Vxx + ^, 



3u 



V'XX 



4w' 

^32 = Vxx + fll — ; 



Vt 



= f2lU^/^V-^' + /22^-"^ + /23^, 



Ut = (V^^'^U: 



3ai/22^i/2 ^ /11/22 



4/21 /21 

^i = /2lW^'/' + /22^'/' + /23, 



= (w^/2^x):r + /ll^^"^/^^^^^+/l2W-^/2 + /i3W, /iss = Uxx + ^ 



2u' 



/?33 = Vxx + ai— ; 



^34 = Uxx + ^iu + /22//21), 



/?34 = ^^x' 



4:V 



aiv. 



Ut = Uxx + fl{u,v), /l35 = + 'y)+aiW + a2, 

aifiuiu,v)u+a2fiu{u,v)-aifi{u,v) 

, P35 = U; 



Vt 



Ut = Uxx + ufi{v) + aiuln(u), 



Vt = as 



fijv) + ai{ln{u) ^ 1) 
ufiiv) 



faiv + a2\ 



Vt = f2iv + /22e-^23" + /: 



24, 



he = Uxx + + aiuln(u) + 

+a2e''i^u + as, 

1^36 = 0; 



/i37 = ^x.+ain, + /i(^-^-^J+ 



23 



P37 = aiVx+ f2iv + f22e^''''+ 
, , , . , N , /24(a3^+a4) 

+724 +(03^ + 04) ?^ + - 



21 
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Ut = Uxx + /2i(a2 - l)w ln(w) + ufi (^-^^ , 

Vt = f2lV + f22U''\ 

hss = + aae^'^^'^^'^^u^ + /2i(a2 - l)uln(u) + + 

= + /i('it) + ^{{n + l)v - ciiu"+'), 
ai 

= ^ii(^X)x + + a2U + a3)/;(w) ^ f (n + 1 - dm)^ + -fiiu)] , 

n + 1 n + l\ ai/ 

/?39 = 0; 

= u^x - - (n + 1)?;) + 

a2\ . / X C^l^S n+l r, n+2 o ^ 

H — uv^ — 2aiv ^ — 2a2V — 03—, 

aiJ n + l u 

h40 = Uxx --^{diuh"^^ - {n+l)v) + fi{u) + {aiv'^^+a2)u+ as, 
diu 

= 0, 

where fij and ak are arbitrary constants. 

Moreover, we found all systems (4.1)-(4.2) which have solutions of the third 
order 

h = Uxxx + hi{t, X, U, V, Ui, Vi,U2, V2), 
P = giVxxx + Plit, X, W, V, Ui, Vi,U2, V2) 
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to the linear determining equations (4.3)-(4.4). The fuU hst includes ten sys- 
tems. We give these systems together with solutions of the linear determining 
equations: 



Ut 

Vt ■■ 


= {uUx)x + hlU^ + fuU + /l3, 
= diVxx + /21W + f22V + /23, 


/i41 


— Uxxx ~l~ fllUx/'^i 
= Vxxx + /ll^ar/2; 


Ut 

Vt ■■ 


= {vUx)x + {hlV + fl2)u + fizV + /i4, 
= diVxx + /21W + f22V + /23, 


h42 

Pa2 


~ Uxxx ~l~ fllUx/'^i 
= Vxxx + /ll^ar/2; 


Ut 

Vt ■■ 


= Uxx + hlU + fl2V + /l3, 

= di{yVx)x + dif2iv'^ + /22?^ + /23W + /24 


h43 


= '^^a;a;a; + /2lWa;/2, 
= 'i^a;a;x + f2lVx/'^\ 


Ut 

Vt ■- 


= Uxx + /uW + /l2?; + /l3, 

= C?i(ui;a;)a; + (/21U + 722)?; + /23W + /24, 


^44 

Pu 


= '^^a;a;a; + /2lWa;/2, 
= 'i^a;a;ar + f2\Vxl'^\ 


Ut 

Vt ■- 


= + f2lU^ + /12W + /13?; + /l4, 

= C?i(ui;a;)a; + {dif2lU + /22)?; + /23W + /24, 


^45 

/545 


= '^^a;a;a; + /2lWa;/2, 
= 'i^a;a;ar + f2\Vxl'^\ 


Ut 

Vt ■- 


= + {f2lV + /i2)w + /13U + /14, 
= C?l(?^?^a;)a; + dif2iv'^ + /22?^ + /23W + /24, 


/l46 
/546 = 


= Wa;a;a; + /2lWa;/2, 
= ^a;a;ar + /21?^a;/2; 


ut-- 

Vt = 


= {uUx)x + f2lU^ + h2U + h2V + fiA, 

= di {yVx)x + dif2iv'^ + f22V + /23W + /24, 


/i47 

;547^ 


= Wa;a;a; + /2lWa;/2, 
= ^^a;a;x + /21?^a;/2; 


ut-- 

Vt = 


= ivUx)x + if2lU + fl2)v + fisU + fu, 

= di{uVx)x + idlf2lU + f22)v + f23U + f24, 


^48 
/548 


= Wa;a;a; + /2lWa;/2, 
= 'i^a;a;ar + f2\Vxl'^\ 


Ut 

Vt ■- 


= {uux)x + hiu^ + /12W + h{v), 

= U21U + h2h{v) + hz)/f[{v), 


^49 — Uxxx + /llWa;/2, 

. ^ ./^'(^) 
/l(^) 


Ut 

vt-- 


= {vU^)x+{fllV + h2)u + fuV + /i4 , 
= /2lW + /22V + /23, 
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^50 


— Uxxx ~l~ /l 1 ^.T / 2 , 
= ^^ara;x + /ll^a;/2. 



As some illustration of the basic method of computing solutions, we consider 
the system 

Ut = {uux)x + — 3u-\-v, (4.5) 
vt = ivvx)x + 2i;^ -2v + u, (4.6) 
with differential constraints 

'^xxx ~l~ '^x ~ 0? '^xxx ~l~ "^a; ~ 0-5 (^•'^) 

It follows from (4.7) that the functions u and v have the representation 

u = wi(t) sin(a;) + U2{t) cos(x) + U3(t), v = 'yi(t) sin(a;) + V2{t) cos(x) + v-^it). 

Substituting this representation into the system (4.5)-(4.6), we obtain ordinary 
differential equations 

u\ = 3ui(w3 -1)+Vi, 

u'2 = 3^2(^3 - 1) + '^2, 
U3 = U3(2W3 - 3) + Ui + U2 + v^, 

v[ = vi{3v3 - 2) + wi, 

V2 = V2iSv3 - 2) + U2, 



5. Invariant solutions under involutive distributions. 

In this section we introduce invariant solutions under involutive distributions. 
Suppose that a collection of p vector fields 

n 

Xs = Y.€s{^)dx, 

i=l 

is given on an open set U C R^. If this collection is linearly disconnected, i.e., 
the rank of the matrix |i^s(a;)| equals p for all x ^ U and satisfies the involution 
condition 

[Xi,Xj] = E 4(x)Xfc, V 1 < ij < p, (5.1) 

k=l 
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where cfj are smooth functions, then this collection generates an involutive p- 
dimensional distribution Dp. A collection of vector fields with these properties 
is called an involutive basis or just a basis. It is well known that a distribution 
Dp is involutive if and only if it possesses at least one involutive basis. 

Definition. A solution u = (/9 to a system of partial differential equation E is 
invariant under an involutive distribution Dp if Dp is tangent to the manifold 
S = {{x,u) : u = (p{x)}. Obviously, the invariance of a solution under Dp 
amounts to its invariance under the operators of an arbitrary involutive basis 
for Dp. 

Now, consider the system of evolution equations 

u\ = F'^{t^x^u^Ua)j z=l,...,m, (5-2) 

where t and x = {xi, . . . , Xn) are independent variables, u^, are functions, 
u = {u^, w"^), and Ua stands for various partial derivatives with respect to 
xi, . . . ,Xn- Denote the total derivatives with respect to t and Xi by the symbols 
Di; and D-j^.. 

Let J^{U,R^) be the space of k-jets on U C R^. Recall that a manifold 
H C J'^ (R^'^^ , R^) , defined by the equations 

h^{t,x,u,up) = 0, j = l,...,s, (5.3) 

is an invariant manifold for (5.2) if the following identity holds on the set 

[^]n[^]: 

Dth^ = 0. 

Here [E] and [H] stand for the differential consequences of (5.2) and (5.3) with 
respect to xi,...,a;„. Denote the involutive distribution generated by vector 
fields Xi, ...,Xr by < Xi, ...,Xr > . 

Lemma 1. Suppose that vector fields 

n m . 

Xk = J2Ckit,x,u)d^^ + J2lk(t^x^u)du,, k = l,...,n, (5.4) 

i=l 3=1 

generate an involutive distribution and that det{^l) ^ 0. If the manifold defined 
by the equations 

K = Ea<. -rfk = 0, 1 < i < m, 1 < A; < n, (5.5) 

i=l 
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is invariant with respect to (5.2) then system (5.2) has invariant solutions under 
this involutive distribution. 

Proof. Write down the cohection of the fields Xi, X„ in vector form as 
follows: 

Acting by the matrix on X, we obtain the involutive collection 

Z = dx + f]du, 

where fj = ^"^r/. The distribution < Zi, Z„ > is involutive. 

The invariant solutions under < Xi, X^ > must satisfy (5.5). The invari- 
ant solutions under < , . . . , Z„ > must satisfy the equations 

=fii{t,x,u). (5.6) 

Obviously, (5.5) and (5.6) have the same solutions. Since Z is an involution 
distribution, the Poisson bracket [Zi,Zk\ vanishes. Consequently, we have 

which means that the consistency conditions for (5.6) are satisfied. 

Using (5.6) and inserting the derivatives of the functions with respect to 
Xk in the right-hand side of (5.2), we come to the system 

ul = G\t,x,u), j = l,...,m. (5.7) 

By the Frobenius theorem, the system of (5.6) and (5.7) is compatible if the 
relations 

D,,G^ = Dtf^l j = l,...,m; k = l,...,n (5.8) 

are valid by virtue of (5.7) and (5.8). Validity of these conditions follows from 
the invariance of (5.5) with respect to (5.2). Indeed, this invariance means that 

A (4. - Vi) = D.r.F^ - Dtfli = 0- (5.9) 

Inserting the derivatives with respect to Xk in (5.9), we see that (5.9) coincides 
with (5.8). 

Remark. If an involutive distribution is generated by analytic vector fields 
Xi, ...,Xp, where p < n, (5.2) is a system of first-order equations with analytic 
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right-hand sides, and the rank of the matrix equals p, then (5.2) has an 
invariant solution relative to Xi, ...,Xp. The proof is carried out by the above 
scheme, but instead of the Frobenius theorem we should use the Riquier theorem 
on the existence of analytic solutions to an autonomous system with analytic 
right-hand sides [5] . 

To exemplify the application of a distribution to constructing solutions, con- 
sider the equation 

z., = Aln«. A = —^ + —^, (5.10) 

which arises in various application [13, 14] and possesses an infinite-dimensional 
algebra of point symmetries [15]. Some exact solutions to this equation can be 
found in [16, 17]. We give a solution to this equation which is invariant relative 
to the pair of commuting operators 

Xi = dx — {u^ + {tu^ — xu^ + u) tan(t))(9u, 

X2 = dy — {tu^ + u — xu^)du- 
The corresponding manifold for these vector fields is 

Ux + u^ + {tu^ - xu^ + u) tan(t) = 0, (5.11) 

Uy + tu^ -\-u- xu^ = 0. (5.12) 

It is easy to verify that this is an invariant manifold for (5.10). Note that the 
vector fields Xi and X2 do not belong to the algebra of symmetries of (5.10). 
The general solution to (5. 11), (5. 12) and (5.10) has the form 

^ A[exp{{x — t)tdLn{t) -\- y]cost + X — t' 

To use vector fields and distributions, we need a method for finding them. 
The classical approach to constructing vector fields relative to which the given 
differential equations are invariant was proposed by S. Lie. A modern exposition 
with many examples and new results was given by L. V. Ovsyannikov [6]. 

A determining equation like (2.5) enables us to find differential constraints 
compatible with the original equation. In the case of differential equations 
in more than two independent variables, we can propose systems of defining 
equations which would enable us to find involutive distributions. 
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Consider the system of involution equations (5.2) and the manifold in J^{U, BI^) 
defined by 

M = <+4(^'^'^) = 0' (5-13) 
where i = 1, m, and j = 1, n. 

Theorem. Suppose that the manifold (5.13) is invariant under the system 
(5.2) whose right-hand sides are polynomials in derivatives whose coefficients 
depend OYit,xi, ...,Xn and u^, ...,u^. Then the functions Mj satisfy the following 
system: 

Dth] + mij{h)\[E] = 0, 1 < i < m, 1 < j < n. (5.14) 

Here mij{h) is some operator representing a polynomial in /if, D^^^hf, Dx^^hf, V^hf 
(/c = 1, m, / = 1, n). The operators mij(h) vanish whenever all hf are zero. 

Proof. We first show that the total derivative of /i*- with respect to t is 
representable as 

Dth] = m^j{h) + -fij, (5.15) 

where niij are operators whose shape is described in the theorem and 7jj are 
functions which may depend only on t, x and u. 
The following identities are valid on [E] : 



A/.;. = i^.F^ + ^+EF^^. (5.16) 



Q\s\„.k 



„ 'U 

Let -r-^f — be a derivative of maximal order on the right-hand side of 

OXi ■ ■ ■ oxfi" 

(5.16) and Sp for some p. By (5.16) and the assumptions of the theorem, 
this derivative enters (5.16) polynomially. Using (5.13), we can write down this 
derivative as follows: 

DZ ■ ■ ■ Dll' ■ ■ ■ Dll{hl) - Dl\ . . . Dll' . . . Dll{gl). 

Note that the second summand involves no derivatives of order |s| and is a 
polynomial in derivatives. Thus, all derivatives of maximal order on the right- 
hand side of (5.16) can be expressed in terms of the total derivatives of the 
functions /i^(r = 1, m, and q = 1, n). Afterwards, it is possible to express 
the derivatives of order \s\ — 1, etc. down to the first-order derivatives. 
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We are left with demonstrating that the functions 7,;^- in (5.15) are all zero. 
By the conditions of the theorem, the manifold (5.13) is an invariant manifold 
for (5.2). Consequently, the following identity holds on [E] r\[H]: 

mij{h)+-iij = Dth] = 0. 

Since the rrii/s vanish on [H], the functions 7^^ are zero on [^]n[-f^]. Once 
the rrii/s are independent of the derivatives of the functions u^, all are 
identically zero. 

Remark. As we see from the proof of the theorem, the choice of the operators 
rriij could be not uniquely defined. 

For example, consider the second-order equation in three independent vari- 
ables: 

Ut = G = F\^:, + F\y + F\l + FWy + F^ (5.17) 
where are some functions depending on u. Suppose that 

hi = Ux + giit, X, y, u) = 0, h2 = Uy + g2{t, x, y^u) = (5.18) 

define an invariant manifold for (5.17). To derive a system of determining 
equations like (5.14), we express the derivatives Dthi and Dfh2 in terms of 
hi,Dxhi, Dyhi, Dlhi,DxDyhi, and Dyhi [i = 1,2). By (5.17), the following 
holds: 

Dthi = D^G + ^ + -^G. 

ot ou 

It is easy to verify that the right-hand side of the last equality is representable 
as 

mii(/ii, /12) = Gu^Dlhi + Gu,,Dlhi + [Gu^ + Dx{GuJ]Dxhi + Gu,Dyhi+ 

+D,[GuJDyh2 + [Gu - DliGuJ - DliGuJ + n]h, + s,h2 + 71, (5.19) 

where ri,si, and 71 are functions depending on hi,h2, and G. Since (5.18) is 
an invariant manifold, the function 71 equals 0. Consequently, the first defining 
equation has the form 

Dthi = mii(/ii,/i2). 
To obtain the second defining equation 

Dth2 = mi2(/ii,/i2), 
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we should replace hi in (5.12) with /i2, x with ri with r2, and Si with S2. The 
following lemma asserts that, under some conditions, solutions to equations like 
(3.13) enable us to construct differential constraints compatible with the system 
of evolution equations (5.2). It is worth to note that the form of the operators 
WLij is inessential, provided that only mij{0) = 0. 
Lemma 2. Suppose that the functions 

n 

h) = E Cjit, X, u)ul^ - g]{t, X, u) 

s=l 

satisfy a system like (5.14) on [E] with fnij{0) = 0. If the vector fields 

n m 

generate an involutive distribution and det(^j) then there is a solution to 
the system of (5.2) and the equations 

/ij = 0, z = l,...,m, j = l,...,n. (5.20) 

Proof. Since the functions hj satisfy (5.14), in view of mij{0) = (5.20) defines 
an invariant manifold for (5.2). To complete the proof, it suffices to refer to 
Lemma 1. 

Finding solutions to general nonlinear equations (5.14) might represent a 
very complicated problem. To simplify the problem, we remove all terms non- 
linear in hf from the operators rriij as it was done above in the case of an 
evolution equation with one space variable. In result, we obtain some linear 
equation 

Dth) + kj{h) = 0. 

As we done above, multiply the coefficients of the operators lij by undetermined 
constants and write down the resultant equations as 

Dth] ^ L,j{h) = (5.21) 

calling them linear determining equations (LDEs). For example, the LDEs for 
(5.17) have the form 

Dthi = Lii(/ii, /i2) = aiGu.^Dlhi + a2Gy^^Dy/ii+ 
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+[a^Gu, + a/^D:,{Gu,^)]Dxhi + a^GuyDyhi + aQDx{GujDyh2+ 

+[a,Gu + a,Dl{GuJ + a,Dl{GuJ]hu (5.22) 

where Li2{hi,h2) is obtained from Lii{hi,h2) by replacing /ii with /i2, a: with 
y, and with 6^. 

Although the above arguments were for systems of evolution equations, we 
can try to extend them to a more general situation. Assume given a system 

ni{u) = F\t,x,u,Ua), i=l,...,m, 

where are linear differential operators with constant coefficients and the right- 
hand sides are similar to those in the case of evolution systems (5.2). To find 
the functions hp we suggest using the following equation in place of (5.21): 

Nim) + Lij{h) = 0, (5.23) 

where the operators Ni are obtained from by replacing partial derivatives 
with total derivatives. Alongside (5.23), it is useful to introduce the following 
analog of B-defining equations [18]: 

Niih^) + Lijih) + E #i = 0, (5.24) 

l<l<m 
l<k<n 

where 1 < i < m,l < j < n, and bfj are functions that may depend on t, x, 
and u. 

We call equations of the form (3.23) quasilinear determining equations (QDEs). 
We exhibit an example of QDEs in finding involutive distributions. Consider 
one of the nonlinear dispersion models describing the propagation of long two- 
dimensional waves [22]: 

T]tt = gdAf] + -j^m + -fi'AT]^ 

where T]{t,x,y) is the deviation of a fluid from an equilibrium state, d is the 
depth of an unperturbed fluid, and g is the free fall acceleration. By translations 
and dilations, we can reduce this equation to the form 

utt - A(uu) - uAu - {Vuf = 0. (5.25) 
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In accordance with the above method, the QLEs for (5.25) have the form 

D^hi - D^Dlhi - D^Dlhi + aiu[Dlhi + D^/ii) + a2UxDxhi + a^UyDyhi 

+a4UxDyh2 + {a^Au + uqUxx + aiUyy + n)hi + gi/i2 = 0, (5.26) 
D^^h2 - DlDlh2 - DlDlh2 + biu{Dlh2 + Dlh2) + b2UyDyh2 + huxDxh2 

+b4UyDxhi + (65AU + 66Wa;a; + ^T^yy + r2)/i2 + ^2^^! = 0, (5-27) 

where and hi are constants, and and qj are functions which may depend on 
t, X, y, and u and which should be found together with hi and h2. The scheme 
for solving (5.26) and (5.27) is completely analogous to the standard scheme of 
the group analysis of differential equations [6, 19]. For this reason, we omit all 
intermediate computations and set forth only the final results. 

If hi and h2 are sought in the form corresponding to the point symmetries 

hi = Clut + ^lux + ^luy + T]^, 

h2 = ^fut + ^2Ux + ^iuy + r/^ 

where and r/-^ are functions of t, x, y, and u, then under the condition (^j)^ + 
(^3)^ + (^1)^ + (^2)^ 7^ equations (5.26) and (5.27) can be shown to have 
solutions leading only to admissible operators for (5.25). There appear new 
solutions only when 

hi = Ux + gi{t, X, y, u), h2 = Uy + g2{t, x, y, u). 

The final form of gi and ^2 is as follows: 

gi = Six + S2y + S3, g2 = S2X + S4y + ss- 

Moreover, the functions Sj (z = 1, 5) depend only on t and satisfy the follow- 
ing system of five second-order differential equations: 

"2 2 
Si + 3si + S1S4 + 2s2 = 0, 

s'2 + 3siS2 + SS2S4 = 0, 
S3 + 3siS3 + 2S2S5 + S3S4 = 0, 
S4 + S1S4 + 2sl + Ssl = 0, 
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^5 + SlS5 + 2S2S3 + 3S4S5 = 0. 

For completeness of exposition, we write down the constants ai and hi {i = 
1, 7) and the functions rj and qj (j = 1, 2) in (3.26) and (3.27) corresponding 
to gi and g2'. 

ai = 61 = 04 = 64 = —1, a2 = 62 = ^3 = &3 = —3, 

= aQ = a7 = h = be = bj = 0, 
ri = 3si + S4, r2 = si + 3s4, gi = 2si, ^2 = 2s2. 
The functions hi and /12 generate the differential constraints 

Ux + SiX + S2?/ + S3 = 0, 

+ S2X + S4?/ + S5 = 0. 

These constraints enable us to find the following representation for a solution 
to (5.25): 

u = — S2xy S3X - s^y + sq. 

Inserting this in (5.25), we obtain the following equation for sq: 

Sq = 3Si + 2siS4 — SiSq + + S3 + 3S4 — S4S6 + S5. 

The system of the six differential equations in the six functions Si deserves 
further study. For example, it would be interesting to find a solution expressible 
via elementary functions. 
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